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ABSTRACT: 

Quantum decoherence is the loss of a system’s purity due to its interaction with the 
surrounding environment. Via the AdS/CFT correspondence, we study how a system deco- 
heres when its environment is a strongly-coupled theory. In the Feynman-Vernon formalism, 
we compute the influence functional holographically by relating it to the generating func- 
tion of Schwinger-Keldysh propagators and thereby obtain the dynamics of the system’s 
density matrix. 

We present two exactly solvable examples: (1) a straight string in a BTZ black hole 
and (2) a scalar probe in AdS;. We prepare an initial state that mimics Schrédinger’s cat 
and identify different stages of its decoherence process using the time-scaling behaviors of 
Rényi entropy. We also relate decoherence to local quantum quenches, and by comparing 
the time evolution behaviors of the Wigner function and Rényi entropy we demonstrate 
that the relaxation of local quantum excitations leads to the collapse of its wave-function. 


Contents 


fo) 


Q W è 


Motivation and Summary 


Formalism for Quantum Decoherence and its Holographic Version 
2.1 Quantum Decoherence 
2.2 Feynman-Vernon and Schwinger-Keldysh 


2.3 Holographic influence functional 


Propagating function 
3.1 Langevin equation from propagating function 


3.2 Evaluating the propagating function and master equation 


Characterization of Quantum Decoherence 

4.1 Preparation of initial state 

4.2 Wigner function: quantum distribution function of phase-space 
4.3 Purity/Rényi entropy 

4.4 Relation to local quantum quench 


Holographic Quantum Decoherence: Case Studies 
5.1 Two Cases 
5.1.1 Straight string in BTZ black hole 
5.1.2 Scalar probe in AdS; spacetime 
5.2 Results and remarks 
5.2.1 String probe 
5.2.2 Scalar probe 
5.3 Initial state profile dependence of decoherence time 


5.4 Numerical studies for discrete time systems 
Conclusion 

Evaluating the propagating function: a review 
Derivation of the master equation for Îsys(t) 


Normalization of fsys 
C.1 Evaluating the normalization A(t) using Trfsys(t) = 1 


C.2 A consistency check: Evaluating the normalization through functional deter- 


minant 
C.3 Another consistency check: The behavior of J(t) and fgys(t) as t > 0 


41 
43 


1 Motivation and Summary 


Quantum decoherence is defined as the loss of ‘coherence’ of a quantum state, and its 
resulting transition to a classical state. What is its mechanism? Is it instantaneous or 
gradual? If latter is it possible to reduce its speed? The development of the quantum 
information science during the last twenty years has not only lifted these questions from 
the philosophical realm to the physical one but also made them pressing issues since a 
realistic quantum computer requires the qubits to remain coherent long enough for their 
operations to complete.! 

The Environment-Induced Decoherence developed by Zurek |1] differs from earlier 
mechanisms of decoherence such as ‘Copenhagen Interpretation’ [2] and ‘Many Worlds’ 
[3] in that it is defined purely from quantum mechanics itself: the notion of classicality, the 
boundary between classical and quantum, the process of decoherence can all be defined and 
computed within the framework of Quantum Mechanics. Adding to its theoretical elegance 
are numerous experimental supports such as [4]. 

In the framework of Environment-Induced Decoherence, the quantum decoherence 
problem is characterized by three elements: the ‘system’ of interest, its ‘environment’, 
and the interaction between the two. Both the ‘system’ and ‘environment’ are quantum, 
together they form a closed quantum system which evolves unitarily. However, the ‘sys- 
tem’ by itself is open: it starts to decohere the moment the interaction is turned on and 
eventually loses all its coherence and becomes a classical state. The ‘environment’ plays 
two crucial roles. First, it determines which states the ‘system’ can decohere into. Second, 
it causes the decoherence to happen via its interaction with the ‘system’. 

The natural question is then how to compute the decoherence rate for different com- 
binations of ‘system’, ‘environment’, and the interaction between the two.? The canonical 
formalism to study decoherence in Environment-Induced Decoherence scheme is Feynman- 
Vernon path-integral [5, 6]. Given the difficulty of decoherence problem (a non-equilibrium 
process of an open ‘system’ possibly at finite temperature), most of the early studies in this 
formalism is limited to the toy model in which both the ‘system’ and the ‘environment’ are 
simple harmonic oscillators (SHO), for which the decoherence rate can be computed exactly 
I7, 8].3 

However, for non-trivial combinations of ‘system’ plus ‘environment’ with slow deco- 
herence rates we need to extend our search beyond this toy model. The Feynman-Vernon 
formalism is a first-principle method that requires no assumption on the detailed nature of 
the ‘system’ or the ‘environment’, therefore is capable of dealing with generic non-trivial ‘en- 
vironment’. This leads to the focus of the present paper: Using Feynman-Vernon formalism, 
we will study quantum decoherence when the ‘environment’ is a non-trivial strongly-coupled 


1 Quantum decoherence is also crucial in resolving the information paradox of black holes and studying 
the primordial cosmic fluctuations during the inflation era of the Universe. 

In general, to mitigate the effect of decoherence on quantum computation one needs to implement 
fault-tolerant algorithms; however here we are only interested in the physical aspects of decoherence. 

3There have also been various approximations and reductions starting from the Feynman-Vernon for- 
malism in order to treat more non-trivial ‘environments’ but they are more ad hoc and their validity need 
to be examined case-by-case (see e.g. [9]). 


quantum field theory that has a dual description in terms of a gravity theory. More pre- 
cisely, we will consider conformal field theories that are dual to gravity theories living in 
anti de Sitter space. In this Gauge/Gravity (or AdS/CFT) correspondence, the strongly- 
coupled field theory can be much more easily studied in terms of its weakly-coupled gravity 
dual. 

The motivation is two-fold. One comes from experiments: the qubits (the ‘system’) in 
the lab are usually embedded in a strongly-coupled condensed matter ‘environment’. For 
instance, a very promising realization of quantum computation is the topological quantum 
computation |10, 11], which utilizes the topologically-ordered phases of certain condensed- 
matter systems. Most of the topologically ordered phases |12] are based on the gapped 
states and hence are robust against ordinary environmental disturbance. However, these 
gapped states or the localized zero-modes such as Majorana fermions [13, 14] (viewed as 
‘system’) still interact with the gapless states on the edge and might decohere via this 
interaction; and the gapless state should therefore be viewed as the ‘environment’ and they 
are usually described by strongly-coupled 2D conformal field theories. 

The other motivation is theoretical. Even beyond the class of the topological quantum 
computation, in general we should search among non-trivial (very likely strongly-coupled) 
theories for decoherence-suppressing ‘environments’. However, quantum decoherence with 
non-trivial ‘environment’ is a difficult problem. Therefore, we could start from those with 
weakly-coupled gravity duals in order to gain some insight to this problem. 

In Feynman-Vernon formalism, all effects of the ‘environment’ can be packaged into a 
certain functional of the fields of the ‘system’ called ‘influence functional’. For the holo- 
graphic quantum decoherence, the key observation is that the influence functional is nothing 
but the generating function of non-equilibrium Green’s functions (Schwinger-Keldysh prop- 
agators) in the ‘environment’ (for the case of linear coupling). This generating function can 
then be easily computed from the gravity side via AdS/CFT correspondence. 

This holographic quantum decoherence applies to any holographic ‘environment’. In 
this paper, we demonstrate its power with two cases of ‘system’ plus ‘environment’: (1) the 
holographic dual of a straight string in a BTZ black hole and (2) that of a scalar probe 
in AdSs. In both cases we solve the dynamics of quantum decoherence exactly. Using two 
quantities (the negative part of Wigner function and the second order Rényi entropy) to 
characterize the quantum decoherence, we then describe the full decoherence process and 
further distinguish its different stages according to the scaling behavior of the second order 
Rényi entropy. We also write a python code to numerically study more complex cases. 
These results not only allow us to understand quantum decoherence in more details; more 
importantly, with our scheme one can study holographic ‘environments’ systematically, thus 
provide valuable insight for the construction of robust qubits. 

We also notice the similarity between the environment-induced decoherence and the 
local quantum quench. We find that the decoherence and the relaxation following the 
quench occur around the same time, and match the scaling behaviors (with time) of the 
entanglement entropy in the quantum quench process with that of the second order Rényi 
entropy in decoherence. This suggests that the relaxation following the quantum quench 
occurs when local excitations decohere. 


This paper is organized as follows. In Section.2 we review quantum decoherence and 
Feynman-Vernon method, and explain how to compute the influence functional from the 
gravity side. In Sec.3 we compute the propagating function and use it to derive the master 
equation for the reduced density matrix. (The detailed derivations are contained in the 
Appendices.) In Sec.4 we study Wigner function, the second order Rényi entropy, and 
relation of our setup to the local quantum quench. Sec.5 examines the two exactly solvable 


cases. Sec.6 contains a summary and discussions. 


2 Formalism for Quantum Decoherence and its Holographic Version 


2.1 Quantum Decoherence 


If our world is intrinsically quantum mechanical, why does it appear to be classical most of 
the time? Can the concept ‘classical’ be defined in a purely quantum mechanical framework? 
How does a quantum mechanical system lose its ‘coherence’ and become classical? Is the 
process of ‘decoherence’, i.e. the quantum-to-classical transition, instantaneous or gradual? 
If latter how can we reduce its speed in order to build a real quantum computer? 

These are the questions a framework for quantum decoherence need to address. His- 
torically, the most prominent schemes are the following two (for a review see |1]). 


1. Copenhagen interpretation [2]: Our world is divided into classical and quantum. A 
measuring apparatus is macroscopic and classical. A quantum mechanical system 
loses its coherence (i.e. wave-function collapse) the moment it is probed by a clas- 
sical measuring apparatus. In this interpretation, the classicality and the boundary 
between the classical and the quantum cannot be defined within the framework of 
quantum mechanics but has to be introduced from outside. 


2. Many Worlds [3]: Unlike in Copenhagen interpretation, the world does not have a 
priori classical subsystems and always evolves unitarily. The decoherence is caused 
by selecting a particular observer or subsystem. 


The ‘Environment-Induced Decoherence’ developed by Zurek |1] differs from these ear- 
lier attempts in that it does not need a deus ex machina outside quantum mechanics to 
come in and announce the transition from quantum to classical. All questions asked in the 
beginning of this section can be answered within the framework of quantum mechanics. 

First of all, a closed system is always quantum. Its state evolves unitarily according 
to the Schrödinger equation, and would never collapse into a classical state by itself. Now, 
let us divide this closed system into a subsystem (which we refer to as the ‘system’) and 
its complement ‘environment’, and we allow the interaction (hence the flows of energy and 
information) between the two. Now the ‘system’ is an open system and due to its interaction 
with the ‘environment’, in general it would not evolve unitarily — even though the total 
system (‘system’ plus ‘environment’) still does. This non-unitarity, hence the decoherence 
of the ‘system’, is caused by the (non-unitary) influence of the environment via a ‘leaking’ 
of the quantum information from the ‘system’ into the ‘environment’. 


The ‘environment’ is also quantum mechanical (unlike in Copenhagen Interpretation) 
and plays two roles. First, it selects from the ‘system”s Hilbert space a subspace of states 
which are stable against the disturbance from the ‘environment’. These states are classical 
states (also called pointer states) — this is the definition of classicality. Second, by inter- 
acting with the ‘system’ the ‘environment’ causes the ‘system’ to lose its coherence between 
these pointer states. The process of decoherence is gradual and its speed can be computed 
given the ‘system’, the ‘environment’, and the interaction between the two. Finally we can 
engineer one or all of these three aspects in order to reduce the speed of decoherence — 
thereby improving the robustness (against decoherence) of the quantum computer. 

To discuss the quantum decoherence process, we need a language to describe the quan- 
tum and classical states in a unified way. The density matrix p is such one. In the basis 
of pointer states, the quantum coherence manifests itself in the presence of ’s off-diagonal 
elements, whose disappearance signifies the quantum decoherence process. At the end of 
decoherence, a quantum state becomes a classical one, whose density matrix loses all off- 
diagonal elements. 

Therefore we need to study the evolution dynamics of the density matrix psy, of the 
‘system’. First, let us denote the density matrix of the total system (‘system’ plus ‘envi- 
ronment’) by /tot(t). As the total system is always quantum, its density matrix evolves 
unitarily according to the Hamiltonian H;.¢ of the total system: 


Atot (t) = eT Hro (tt) A oe (tije ott) (2.1) 


where /tot(t;) is the initial density matrix of the total system at t = t; and we set h = 1 
throughout the paper. We assume the factorized initial condition, i.e. at t = t; the ‘system’ 
and the ‘environment’ are unentangled: 


Ptot (ti) = Psys(ti) Q Êenv (ti) . (2.2) 


Once ptor(t) is known, the reduced density matrix of the ‘system’ of interest is given by 


tracing out the degrees of freedom of the ‘environment’: 
Paya (t) = Treny Ptot (t) : (2.3) 


Then we can prepare the initial ‘system’ to be in a pure state, and the quantum 
decoherence process of the ‘system’ is encoded in the way the reduced density matrix 
Psys(t) evolves from that of a pure state into a classical one. Formally, the dynamics of the 
reduced density matrix can be encoded in a Schrédinger-like equation, i.e., 


dpsys(t n 
je — ys fost] (2.4) 
where Hys is the renormalized Hamiltonian of the ‘system’. The --- terms characterize the 


non-unitarity of the influence from the ‘environment’ and is responsible for the quantum 
decoherence. In a general lab experiment, the transition of the ‘system’ from quantum 
to classical happens very fast. (Indeed, in general it appears to happen instantaneously. ) 
Now, in the framework of environment-induced decoherence, we can actually compute this 


decoherence time-scale by a direct study of the time-evolution of fsys(t). In particular we 
can verify that when the ‘environment’ is macroscopic as in a lab experiment, and the 
transition is indeed very fast. 

Most of the computation in this paper is carried out in the path-integral formalism, 
therefore let us now look at the Lagrangian description. The Lagrangian for the total 
system (‘system’ plus ‘environment’) consists of three parts: 


Ld, x] = Lsys[¢] + Lenv[x] + Lint[®, x] - (2.5) 


The first two terms Lgys[¢] and Leny[x] define the ‘system’ and ‘environment’, respectively. 
We use ¢ to denote collectively the degrees of freedom of the ‘system’, and y those in the 
‘environment’. The interaction between the two is given by Lint. In this paper, we consider 


the simple case of linear coupling: 


Lintl?, x] = 9 $ Olx] (2.6) 


where g is the coupling constant and O[x] is a given function of x. 

We also assume that the number of degrees of freedom of the ‘system’ is much smaller 
than that of the ‘environment’, such that the back-reaction of the ‘system’ to the ‘environ- 
ment’ can be ignored during the time scale of quantum decoherence of the ‘system’. This 
is a realistic assumption and has been adopted in the past study of quantum decoherence 
(see e.g. [6, 8, 15]) 

The main difference of the present paper from earlier studies is the following. Previous 
studies of quantum decoherence usually consider the case in which the ‘system’ ¢ is a simple 
harmonic oscillator (SHO) and the ‘environment’ y consists of a collection of independent 
SHOs at thermal equilibrium with temperature T = 1/3. The (linear) interaction between 
the two is dictated by the spectral weight C; of the SHOs in the ‘environment’: 


LEM) = 9 6S Cixi - (2.7) 


This model is quadratic therefore the influence functional and evolution dynamics of sys can 
be obtained exactly in this case. Indeed, this is the simplest model to study non-equilibrium 
processes in the presence of environmental influence: the Feynman-Vernon formalism was 
developed using this model [5]; and later it was used to study the quantum Brownian motion 
in [6] (which derived its Langevin equation after making a suitable choice of the spectral 
weight C;); finally it was also used to study quantum decoherence in |8, 16]. See also [17-20] 
for the similar consideration of the fermionic version. 

However, interesting physics happens when we go beyond the models with simple har- 
monics oscillators — this is what we will do in this paper. Both the ‘system’ and the 
‘environment’ have two different aspects: the theory itself (given by the Lagrangians) 
and its physical state (given by its time-dependent density matrix). Since the ‘environ- 
ment’ is assumed to be in a thermostatic state throughout the decoherence process (i.e. 
Îenv = e7? Henw), we will allow the ‘environment’ Lagrangian to be arbitrary, although later 
in the actual computation we focus on theories with a dual gravity description so that we 


can use the holographic machinery to obtain results that would have been hard to compute 
directly in the field theory side. And we are interested in the difference between differ- 
ent environment Lagrangians. On the other hand, we consider the ‘system’ as a probe to 
study the decoherence effects of different ‘environments’; therefore we will choose a simple 
‘system’: a canonical scalar with Lagrangian 


Laysl6] = -300 - FPE. (2.8) 


The influence functional can be regarded as the probe’s effective action, which is obtained 
after the ‘environmental’ degrees of freedom y are integrated out. 


2.2 Feynman-Vernon and Schwinger-Keldysh 


In the coupled total system given by the Lagrangian (2.5), we are only interested in the 
evolution dynamics of the ‘system’, but not in the detailed dynamics of the ‘environment’. 
Therefore the degrees of freedom of the ‘environment’ should be integrated out and its 
whole influence on the ‘system’ packaged into one or a few quantities. For the generic 
time-dependent problem at hand, this task requires the path-integral formalism devised by 
Feynman and Vernon in [5]. 

In this subsection, we review the Feynman-Vernon (FV) formalism. The main point 
of FV formalism is in rewriting the evolution (2.1) in the path-integral representation, and 
thereby integrating out the ‘environmental’ degrees of freedom y to produce an ‘influence 
functional’ that contains all the ‘environmental’ effect on the ‘system’. As the evolution 
of density matrix involves both forward (et!) and backward (e’“**t') propagators, the 
path integral should be formulated on an ordered closed-time path (Keldysh contour kK), 
i.e., from t = t; to ty and then back [21, 22]. The ordinary causal Green’s function is 
replaced by a 2 x 2 matrix of Green’s functions (a.k.a. Schwinger-Keldysh propagators) to 
account for the two branches in the Keldysh contour. The influence functional is precisely 
the generating function of these Schwinger-Keldysh propagators [23]. 

This formalism, including the connection between Feynman-Vernon and Schwinger- 
Keldysh, is valid for any initial density matrix eny of the ‘environment’. For our present 
case of thermostatic environment with Jeny = e~BHenv the effect of the thermo-average can 
also be represented as a path-integral, but along the imaginary time direction from t; to 
ti — ip. The thermo-Keldysh contour C is along the path 


C: ti === t > tr 10 > t,—10 > ti—iß. (2.9) 


In the original thermo-Keldysh contour, ø is chosen to be 0. Later it was shown in [24] that 
g can actually be chosen arbitrarily since the ‘environment’ is thermostatic. To compare 
with the bulk computation, the symmetric choice o = £ is the most convenient one [25]. 
Therefore in this paper we will use the thermo-Keldysh contour C with o = i as shown in 
Fig. 1. 

Now let us use the method of Feynman-Vernon to compute the reduced density matrix 
fsys at the final time tp. Take an arbitrary element (¢4|Asys(tr)|d_). It is given by the 


201708.00144v1 


chinaXiv 


ti, Xa O+ X+ P+ tp Xr P+ 


ti — io, Z4- 


XB 
t; —if 


Figure 1. The thermo-Keldysh contour C. The two horizontal segments represent the closed-time 


path (i.e. Keldysh contour X), which is the contour for a generic ‘environment’ density matrix Penv- 


When the ‘environment’ is thermal with temperature 4, the average of the ‘environment’ can be 


accounted for by the path-integral along the additional vertical segments. 


total density matrix Prot (ty) via (2.3): 


TA ETE, -| te Pere an (2.10) 


where we use the bar to label the final values (at t = ty) of fields. Since we know that the 
total density matrix ftot evolves by (2.1), after inserting complete sets of fields at t = t; we 
can expand (2.10) into 


Gslpostt)I6-) = | dzdz+dz-dġ+ aĝ- 
Galera ti) Noy x4 ) ` (ġ t , X4 |Psys(t i) ® Penv (ti) |¢ xX ) ` (¢ X jetos =t) 15 _ x2) 
(2.11) 


where we used tilde to label the initial values (at t = t;) of fields. Note that the initial 
density matrix of the ‘environment’ feny(t;) can be arbitrary. 


The integrand of (2.11) consists of three terms. The first and third are forward and 
backward propagators, respectively, which can be rewritten in the path-integral represen- 


tation as 
x pt 
(pa, gale Et 1b, k) =f" Do, f * pyy e Kf clero 
eo ee l (2.12) 
(b_, zlet) Do_ Dy_ eift dt Lle- x-] 
$ š- 


where ¥; = X— = X and we have used the subscript +/— to denote fields living on 
forward/backward time path. 

Now let us extract all the information of the ‘environment’ field x from the r.h.s of 
(2.11). Plugging (2.12) back into (2.11) and recalling that the Lagrangian in (2.12) is given 
by (2.5), we can package all the information of the ‘environment’ field x into an influence 
functional F defined as: 


Fibs =f LE ENEE 


(2.13) 
f Drs f Dy_ aa dt( (Lenv[X+]— Lenv[X ]+Lint [64 X-4 [l= Lint |? X D. 
X+ 


Once the influence functional is known, (2.11) is given by 


(Bless) = | dbxdd_J [5.,6-st7182,6-+4] Grt e 


where J is the propagating function that dictates the evolution of sys from t = t; to ty, 
and it depends on the ‘system’ Lagrangian and ‘environmental’ influence functional via: 


$+ $- 


J |, op st ploy ,¢ i = a Dox 3 Do- eile! ar (Leys [$+] Feyslo—] F (b+, o-] i 
$+ $- 
(2.15) 
Therefore the most crucial task is to compute the influence functional F. 
First, (2.13) can be written concisely as 
Flos, ġ-] = (Tet fk FintlOxl) nv (2.16) 


where the correlator is averaged w.r.t. the initial density matrix of the ‘environment’: 


lo. Jeny = Treny [env -..] (2.17) 


and K denotes the closed-time path (i.e. the two horizontal segments in the contour C of 
Fig. 1.), Ti is the path-ordering operator. For a generic interaction Lint, the influence 
functional F is difficult to compute. However, when the interaction Lint is linear in ¢ as 
given in (2.6), (2.16) becomes 


Fip, p-] = (Tke? POM) ny . (2.18) 


If we consider the ‘system’ field ¢ as the source of the ‘environment’ field y (or more precisely 
of O[y]), the influence functional F with the linear coupling (2.18) is precisely the generating 
function of the real-time Green’s functions at finite temperature (a.k.a. Schwinger-Keldysh 


propagators), which are defined as 


Gss (1, 2) = —i(TkOs(1)Os (2)) env (2.19) 
where 7x is the path-ordering operator along the closed-time path K, s,s’ = +, and O,(1) = 
O[xs(t1, £1)] [23]. Namely 

l = 
daaa i l _ (2.20) 


"3bs(1) 8ps (2) '=9 


Given the influence functional, the Schwinger-Keldysh Green’s functions can be ob- 
tained via (2.20). However, for the decoherence problem at hand, the question is the 
inverse: How to compute the influence functional F once we know the Schwinger-Keldysh 
Green’s function Gss (1,2) of the ‘environment’? This question can be answered when the 
coupling between ‘system’ and ‘environment’ is weak (i.e. g < 1), where we can approxi- 
mate by keeping only the quadratic coupling term (in the exponent of F) and dropping all 


higher order terms:* 


2 
—ilnF[é1,¢_]& -E fats f atta! S sgn(ss')ġs(£)Gss (a —2')bg(x'). (2.21) 
s,s! 
In fact, when the ‘environment’ is composed of a collection of simple harmonic oscillators, 
the coupling is necessarily quadratic (in ¢), and (2.21) is actually the exact answer. 

To summarize, the strategy of studying the decoherence (in the FV formalism) is to 
first solve the Schwinger-Keldysh Green’s function of the ‘environment’, then using (2.21) 
we can immediately obtain the influence functional F, which in turn gives the propagating 
function J via (2.15), and then it is straightforward to study the evolution of sys using 
(2.14). 

Before we proceed, let us first simplify (2.21) further. First of all, out of the four 
Green’s functions Gss only three are independent, since by definition 


G+ +G =G +G ` (2.22) 


The remaining three linear combinations of Gss contain two types of information on the ‘en- 
vironment’. The advanced and retarded Green’s functions {Ga, Gr} encode the dispersion 
relations of the ‘environment’, whereas the symmetric Green’s function Gym describes the 
actual state of the ‘environment’. Therefore it is more convenient to switch to a basis for 
the matrix Green’s function such that its four components are manifestly {Ga, GR, Gsym, 0}. 
Accordingly, instead of {¢+4, ¢_}, we should switch to the basis {%, A} in ‘system’ [22, 27]: 


OE T= i a OEF TERE (2.25) 


The field A is the difference (i.e. relative) between the two fields ¢, and ¢_ (from the two 
branches of the contour), and © is the average — therefore {X}, A} is called the ‘ra’ basis, 


yu 


in which the influence functional can be written as 


i 


. ts ts t / 1 1 1 1 
-imf =- f ar f i [ats dix a aaa 
2.26 


‘It was explicitly shown in [26] that the leading term of the influence functional yields the generating 
function for the Schwinger-Keldysh Green’s functions of O[,x] in a thermal reservoir. 
5We remind that {Ga, GR, Gsym} are defined as follows 


: À 1 
iGa = O(t2 = t1) (0 (2), O) env; IGr = O(t1 — t2)((00), O(2))env, Gsym = 5 A001), O(2)})env - 
(2.23) 
and they are related to the Schwinger-Keldysh Green’s functions via: 
Ga = G++- G-+, Gr=G44—-Gy-, Gsym= 5 (G44 Gi) x (2.24) 


—10- 


where we have used the relation 
Ga(1,2) = GR(2, 1). (2.27) 


So far, the derivation applies to generic feny(t;). Indeed, as long as we can obtain 
the ‘environment’s Green’s function, we can study the decoherence in this ‘environment’ 
using (2.26) no matter how exotic the ‘environment’ is — even if the ‘environment’ itself 
is in non-equilibrium (with interesting examples being de Sitter space and non-thermal 
‘environment’). 

Now let’s restrict to the most realistic and simplest situation when the ‘environment’ 
is in a thermo-equilibrium with temperature zi i.e. the density matrix peny = e~ bhen, 
In such cases, the trace in (2.17) reduces to a thermal average and can be accounted for 
by a path-integral from t; to t; — iĝ. Namely, the Keldysh contour K (the two horizontal 
segments in Fig. 1) is supplemented by the vertical segment to give the full thermo-Keldysh 
contour C defined in (2.9) and represented in Fig. 1. The fields that appear in the definition 


of Schwinger-Keldysh propagators now become O+(1) = O[x+ (tı, #1)] with 


x+(t) 2 etHenvt e—tHenvt , x (t) = et Henv(t—t0) . 9—tHenv(t— ia) f (2.28) 


When the ‘environment’ is in a thermo-equilibrium,® its Green’s functions satisfies 
another constraint, namely the Kubo-Martin-Schwinger (KMS) condition [30]: 


G4_(t — ib, Z) = G_+(, 2) (2.29) 
which gives rise to the condition 
Gsym(w) = —[1 + 2n(w)]ImGR (w) (2.30) 


in the frequency domain, where n(w) = oy is the thermal distribution of the ‘environ- 
ment’. Eq. (2.30) relates the symmetric Green’s function to the imaginary part of the 
retarded Green’s function, which is the spectral function for the O[x] excitations. 

To summarize, when the ‘environment’ is thermostatic, all its dynamical information 
can be encoded in the retarded Green’s function GR. Note that there is no quadratic term 
for X in (2.26) due to the condition (2.27). Moreover, the quadratic form of (2.26) ensures 
that one can further solve the propagating function for fgys, i.e. (2.15) in a closed form as 
long as Ssys[¢] is also quadratic. This implies that the dynamics of psy, can be accurately 
determined once Gp is known. 

We consider the relation (2.26) (or together with its higher order term) the most es- 
sential element in the study of quantum decoherence as it connects the Feynman-Vernon 
(the theoretical framework of environment-induced decoherence) to Schwinger-Keldysh (the 
machinery that allows us to capture and compute the influence from the ‘environment’). Its 
interpretation is the following. Once the ‘system’ starts to interact with the ‘environment’, 


® Otherwise, the condition (2.29) and (2.30) are not valid. However, in this paper we will restrict to the 
thermostatic environment at all time as the probe is small and local compared to the environment. For 
discussions of quantum decoherence in non-thermostatic environments, see [17, 18, 20]. 
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it is under the influence of the ‘environmental noise’. The basic properties of this noise, 
such as its energy and lifetime, is characterized by the ‘environment”s retarded Green’s 
function Gr (in terms of its poles and zeros), from which one can extract the corresponding 
transport coefficients of the thermal reservoir by, for instance, the Kubo formula. (This step 
is the standard application of Schwinger-Keldysh formalism.) Then from the viewpoint of 
the ‘system’, the transport phenomena carries away its quantum information and causes its 
(gradual) decoherence. 

At this point, one might wonder why this natural connection hasn’t been very visible 
in the literature of quantum decoherence. There are two reasons. In the usual discussion of 
quantum decoherence, the ‘environment’ is modeled by a collection of SHOs for simplicity,’ 
and its spectral function (i.e. the spectral weight C; in (2.7)) can then be chosen at will 
and put in by hand (instead of being determined from first principle QFT computation). 
Furthermore, the simplicity of the SHO model makes the machinery like (2.26) unnecessary 
since in this case the influence functional can be easily obtained in closed-form. It is 
this simplicity that has delayed the application of Feynman-Vernon formalism to more 
general ‘environments’ and therefore the study of quantum decoherence in more physical 


and interesting situations. 


2.3 Holographic influence functional 


The relation (2.26) is central to Environment-induced decoherence and applies to any ‘en- 
vironment’; now let us apply it to cases beyond the toy model with SHO ‘environment’. A 
crucial step is to evaluate the retarded Green’s function Gp of O[x] (the operator in the 
‘environment’ that couples to the ‘system’ field ¢). After obtaining the retarded Green’s 
function, we can plug it into (2.26) to calculate the influence functional from which one can 
explore the behaviors of quantum decoherence. The retarded Green’s function is in general 
difficult to evaluate except for the free theory. This is part of the reason that the quantum 
decoherence for the interacting ‘environment’ is much less studied. 

However, besides the free theory, there exists a class of strongly-coupled theories whose 
retarded Green’s functions can be obtained in practice. These are the theories with dual 
gravity descriptions. The most developed class of Gauge/Gravity duality is the AdS/CFT 
correspondence, in which the (non-gravitational) field theory is a d-dimensional conformal 
field theory and its gravity dual lives in a (d+1)-dimensional asymptotically AdS space. The 
zero-temperature (Euclidean) two-point functions of these CFTs have power-law behaviors; 
whereas the finite-temperature real-time Green’s functions cannot be determined by the 
conformal symmetry alone, but can be evaluated through its gravity dual. 

In the present context, our probe field ¢ plays the role of the external source coupled 
to the operator O[x]. In the gravity dual, ¢ is the boundary value of a massive elementary 
field in the asymptotically AdS space, with its mass determined by the conformal dimension 


TA notable exception is the spin-bath model studied in [29]. Modeling localized modes such as nuclear 
spins and defects in the environments, a spin-bath is intrinsically strongly-coupled to the ‘system’ therefore 
cannot be mapped to an oscillator bath: whereas the back-reaction from the ‘system’ to the oscillator-bath 
is usually negligible, the coupling between the ‘system’ and the spin-bath is so strong that the dynamics of 
spin-bath is slaved by the ‘system’ [29]. It would be interesting to holographically model the decoherence 
due to a spin-bath. 
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of operator O[x]. Then, the generating function for the correlation functions of operator 
O|x], i.e., the influence functional, is determined by the on-shell action of the corresponding 


massive field [31]: 
FY dg] = (e9 Se ditt @ Obl) — ¢§Sburx[®|¢] (2.31) 


where ®|g is the on-shell bulk scalar with ¢ being its value at AdS boundary, and Spur is 
the bulk action of 6. The superscript A denotes the “holographic". 

The prescription of calculating the on-shell bulk action to derive the Euclidean two- 
point functions was proposed in [31]: one imposes Dirichlet boundary condition at the AdS 
boundary and needs to choose appropriate boundary action when evaluating the on-shell 
action. For evaluating the Lorentzian correlation functions such as the retarded Green’s 
function, especially at finite temperature (corresponding to a black hole in asymptotically 
AdS space), the choice of the boundary condition at the black hole horizon and the corre- 
sponding boundary action are more subtle. The subtlety arises because there is a coordinate 
singularity at the horizon, thus one should choose the appropriate boundary condition for 
| at horizon to have a smooth solution. This choice corresponds to the thermal vac- 
uum for the dual CFT. To compute the retarded Green’s function, one should choose the 
in-falling condition for the on-shell @|4 at the black hole horizon, and the result is [28, 32]: 

I(z,r 
Gi? (a) = = lim n (2.32) 
where zx is the transverse coordinate of the AdS space and r is the radial one with the AdS 
boundary at r — oo. Here H|gjg is the conjugate momentum of ® with respect to the 
on-shell bulk action. 

In fact, the denominator of (2.32) is the leading term in the large r expansion of ®|¢, 
and the numerator is the sub-leading one. After scaling away the r factors in both the 
numerator and denominator of (2.32), it can be seen as a linear response theory for the 
dual CFT: Gy = Ole if one identifies H|g)4 as (O)|¢ after some appropriate rescaling.® 

It was shown in [25] that the on-shell action for a bulk scalar in the AdS black hole 
background yields the same form as the r.h.s. of (2.21). The forward and backward paths 
in the Schwinger-Keldysh formalism correspond to the time path in a causal patch and the 
one in its mirror image of the Penrose diagram for an eternal AdS black hole. 9 Combining 
this result with (2.31), we summarize that the holographic influence functional is given as 
follows: 


f tf ty h a 
—iln F™ = —¢? L dr i dr’ l d'r dîr’ [AGR (w—2')B(2")— A(x) GQ), (e—2')A(2') 
(2.33) 
where Gh is evaluated as in (2.32), and Ge is related to gh by the KMS condition (2.30). 


This is the holographic version of (2.26). Again, the quadratic form of (2.33) is guaranteed 
by the weak coupling g so that the higher order correlation terms are suppressed. 


8For a massive canonical bulk scalar ® in a (d+ 1)-dimensional AdS space, $() = lim, ®¢(x, r) and 
(O)|¢ = lim, r*~“H(a,r)|e\g with A the conformal dimension of O. 

°This kind of identification was first proposed in [33] for Schwarzschild black hole, and later generalized 
to AdS one in [34]. 
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Here, two remarks are in order: 


1. Although one can obtain the holographic influence functional by treating the ‘envi- 


ronment’ as a holographic CFT, one cannot naturally incorporate the kinetic term 
for the probe field in this context as it is treated as an external source. Thus, we need 
to put in probe’s kinetic term such as (2.8) by hand. It would be more satisfactory if 
the dynamics of the probe field could also be embedded in the bulk holographically. 


. The holographic CFTs are usually gauge theories therefore the operator fx] should 


be a gauge invariant observable. Hence the probe can only decohere through a partic- 
ular channel of physical observable to which it couples. The probe and its decoherence 
behavior are characterized by the spectral density of O[x] encoded in the retarded 
Green’s function. This is in contrast to the case of SHO ‘environment’ in which the 
probe SHO couples directly to all the ‘environmental’ degrees of freedom and the 
probe and its decoherence behavior are characterized by the spectral weight C; in 


(2.7). This difference is illustrated in Fig. 2. 


Figure 2. Different decoherence mechanisms for the probe ¢. Left: ¢ in thermal reservoir of simple 
harmonic oscillators. Right: ¢ in holographic thermal reservoir. 


3 Propagating function 


As shown in the previous section, the evolution of the ‘system’s reduced density matrix 


Psys is given by an integral transform (2.14), whose kernel is the propagating function J 


given by (2.15). The propagating function J is determined by two inputs via the path- 


integral (2.15): (1) the data of the ‘system’ (i.e. its Lagrangian and its initial and final 


states {¢ġ 


$ 


}); (2) the effect from the ‘environment’ (i.e. the influence functional F, 


which can be expressed in terms of Schwinger-Keldysh propagators of the ‘environment’). 


In this section, we first combine the two inputs to rewrite the propagating function directly 


in terms of the ‘environment”s Green’s function plus the ‘system”s data. This then allows 


us to derive the master equation of Îsys. 
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3.1 Langevin equation from propagating function 


The propagating function J given in (2.15) can also be interpreted as the path-integral of 
the ‘system’ on the Keldysh contour K: 


J [8+ 9-st/1d4, -iti = A Do4 ia Do elSenl?+9-] (3.1) 


in which the action of the ‘system’ is modified into Seg due to the influence of the ‘environ- 
ment’: 


aad =f” dt (Leysld4] — Loyoló-]) — in Flor, 9-1. (3.2) 


Now we switch to the ‘ra’ basis which is more suitable when studying the non-equilibrium 
process in terms of ‘environment”s Green’s functions, and set t; = 0 and ts = t, with which 
J becomes a partition function of {X, A}: 


5 A 
J [£, A; t£, A; 0] a. pe f DA eiSeal=.Al (3.3) 
x A 
with action 


Sue, A] = [ arb )A(r) — Q E(T)A(T)] 
(3.4) 


=y aK ar f ar dr'[A(T)GR(T — TEI) — LA(t)Geym(t — TACT) - 


where we have used the ‘system’ Lagrangian (2.8) and the expression of F in terms of 
Schwinger-Keldysh propagators (2.26), and = ce , and we suppress the spatial dependence 
for simplicity. 

Before we proceed to evaluate the propagating function J using (3.3) and (3.4), let’s 
first check that it has the correct semi-classical limit. In the semi-classical approximation, 
J should give the Langevin equation for the quantum Brownian motion as shown earlier in 
[6, 35, 36]. Now we review the derivation. Viewing J as the partition function of {%, A}, 
and recall that the ‘average’ field © is slow and heavy whereas the ‘relative’ field A is fast 
and light, we can integrate out the light field A and obtain the classical trajectory of the 
heavy field X. 


First, via a Legendre transform, the A? term in the r.h.s. of (3.4) can be rewritten as 
eae i dr i dr'A(T)Gsym(T=T')A(T') = fo a R drA(r)&(r)— zz Jo dr i dr'€(r )Gaym(r- TJET’) . 


(3.5) 
Namely the fast field A can be regarded as the response field of a Gaussian random force 
€ which has the correlation 


(ERIE? Vaux = g Gsym(T = T’) . (3.6) 


Intuitively, integrating out the ‘relative’ (light and fast) mode A should provide random 
kicks to the ‘center-of-mass’ (heavy and slow) mode £. Therefore we plug (3.5) into (3.3) 
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and (3.4) and integrate out A to obtain 


J= [rete dr f° dr'€(T)Gsym(T—-1'JE(7’) [x ee) | dr!Gp(t—1')=(r')—-€(7)] 
(3.7) 
The argument of the delta-function gives the classical equation obeyed by the center-of-mass 


field X: z 
X(T) +0?d(r) + al dr'Gr(r — 7')X(r’) = E(r) (3.8) 


which is precisely the Langevin equation, where © is the slow field and €(r) the noise that 
satisfies (3.6). The fluctuation-dissipation theorem for this quantum Brownian motion is 
nothing but the KMS condition (2.30). 


3.2 Evaluating the propagating function and master equation 


Unlike the discussion in the previous subsection, to fully characterize the dynamics of sys 
and examine the quantum decoherence of the probe, one should carry out explicitly the 
Gaussian path integral of (3.3). Although the calculation has been done, for example in 
[8], we review it in Appendix A for completeness and to fix the notation. Besides, there 
are some issues about the normalization of sys which were not considered explicitly in [8], 
thus we carry out the explicit calculation in Appendix C. 

Since in the decoherence process, we are mainly interested in the time-dependence 
of Psys, from now on we will focus on the case in which the ‘system’ field has no spatial 
dependence. For example, this is the case when the ‘system’ is a particle and the ¢ field its 
coordinate, or if we only consider one particular momentum mode of ¢. The result of J is 


then summarized here: 


ERISA l i Weal, 5, Š, Ast 
JÈ, A, t|5, A, 0] = TRI Wer[, ee Ae (3.9) 
where 
= s e š So r S $ zÑ s 1 x92 Que =) 
Wes = (HOË + AOE) A- (AOE + fi()E) A-5 (aooA? + (a01 + a1)AA + 011A?) 
(3.10) 
and 


a(t) = èf ar f dT'gilT)Gsym(T — 7’)g;(7’) . (3.11) 


For simplicity, we set t; = 0 and tr = t from now on. The functions {h, fo, f1, go, gı} are 
defined as follows. First, h(7) is the central one based on which all others are defined; and 


it is the solution of the initial value problem: 


A(r) + hr) + 9? 1 dr'GrRr(T—T')h(r)=0 with A(0)=0, A(O)=1 (3.12) 
0 


which can be solved via Laplace transform:!? 


1 
8 +02 + PGp(s) 


10n practice, this is not necessarily easier than directly solving (3.12) numerically unless when G(r) 
can be obtained analytically. 


(3.13) 


h(r) = £7} [h(s)| with f(s) 
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Then { fo, f1, 90,91} are defined in terms of h via: 


AGS -ihr + h(n) flr) = gyn (3.14) 
and 
go(T) = fit-7), g(t) = fo(t—T). (3.15) 


Note that { fo,1, 90,1} are actually functions of both 7 and t, with t entering as the boundary 
point. For details of the derivation above see Appendix A. 

With the above explicit form of the propagating function, we can now derive fgys’s 
master equation (the evolution equation of a density matrix) starting from (2.14). In the 
‘ra’ basis, the l.h.s. of (2.14) is a function of {©, A}, and we compute its time-derivation. 
Since the derivation is rather long, we leave the details in Appendix B and here only write 
down the final result: 


Opeys(S,A,t) . 8 25 nan ex 
1 Ot = | PIN + Qen lt)2A]psys(%, A, t) 
ð : er a 
tin Aa, + vA + i73(t)A*] psys(%, A, t) 
2 
+iya(t) Bea Psy», A, t) (3.16) 


(for detailed expression of Oren and y; please refer to Appendix B). 

Let us now explain the master equation (3.16). In the first line, Qten(t) is the renor- 
malized frequency, and this line is the renormalized Liouville term which preserves the 
unitarity. The terms in the second line are responsible for the quantum decoherence. 7; (t) 
with i = 1,...,4 are functions of time in generic situations; in the special case where all 
y's are time-independent, the evolution dynamics is Markovian. Note that the last term is 
absent in [8], however it can be re-absorbed into the kinetic term of the Hamiltonian via 
a renormalization of the mass. The difference could be due to the different routes in the 
derivation of (3.16): We obtained (3.16) by directly taking time derivative on (3.9), whereas 
in [8] the counterpart of (3.16) was derived via a perturbation analysis of the path integral 
expression of (2.15). 


4 Characterization of Quantum Decoherence 


After deriving the explicit form of the master equation (3.16), we could now solve it to 
obtain the density matrix Îsys(t) for a given initial state. Alternatively, since we have the 
explicit form of the propagating function given by (3.9), we can also directly compute Pgys(t) 
via (2.14). Once we obtain fsys(t), how can we use it to characterize a decoherence process? 

Recall that the ‘environment’ is not only responsible for destroying the coherence in the 
‘system’: first and foremost it selects from the Hilbert space of the ‘system’ a subspace of 
states that are most stable against the ‘environmental’ disturbance; these states are called 
pointer states, and become classical states at the classical limit A —> 0 [1]. 
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In the basis of these pointer states, the quantum coherence is encoded in the off-diagonal 
elements of the density matrix of the ‘system’, and the Environment-Induced decoherence 
is then the process of the gradual decrease (to zero) of these off-diagonal elements. 

However, for our present ‘system’ of a scalar ¢, it is not feasible to directly look at 
the explicit elements of sys, for the following two reasons. First, since @ takes continuous 
values, sys is an infinite-dimensional matrix. Second, “decoherence is the disappearance of 
Psys’s off-diagonal elements" is a basis-dependent statement; however, the pointer states are 
not exactly stable under the ‘environmental’ disturbance (since we do not take h — 0 limit) 
therefore the pointer basis of the density matrix is not invariant during the decoherence 
process. 

To solve the first problem, we can coarse-grain the ‘system’ — i.e. to prepare its initial 
state to resemble a simple finite-dimensional system, even a qubit (the quintessence of 
Schrédinger’s cat). In this section, we will first explain how to prepare (an approximation 
of) Schrédinger’s cat with the scalar ¢ in our ‘system’. 

The second problem can be solved in two ways. The first basis-independent information 
on the degree of ‘coherence’ can be found in the Wigner function W (X, p, t) (the Fourier 
transform w.r.t. the fast variable A = + — ¢_) of a density matrix p(X, A,t). W(z, p, t) 
is the quantum version of the distribution function f(x,p,t) in the phase space {z, p}; 
and unlike f(z, p,t), W(x,p,t) is not always positive-definite, but will become so once the 
system completely decoheres and becomes classical. Therefore, we can unambiguously use 
the disappearance of W (r, p,t)’s negative parts to characterize the decoherence. 

Another way to solve the basis-dependence problem is to look at the entanglement 
entropy or the Rényi entropy since they are scalars therefore are independent of the basis 
choice. In this paper, we will choose the second order Rényi entropy (which is the logarithm 
of the ‘purity’). As we will see, the Wigner function and the second order Rényi entropy 
are both very effective in characterizing decoherence. 

Finally, we will show in this section that our setup to study quantum decoherence is 
almost parallel with the one for the local quantum quench [37|. Thus, it is interesting to 
establish the connection (if there exists any) between quantum decoherence and the local 
quantum quench. 


4.1 Preparation of initial state 


A prototype for the study of quantum decoherence is Schrédinger’s cat. The initial pure 
state of the cat is in a superposition of the two pointer states |Dead) and |Alive):!! |Cat) = 
a (|Dead) + |Alive)). After decoherence, the cat has a classical probability of being 50% 
|Dead) and 50%|Alive) (after she comes into contact with the environment for long enough 
but before we open the window to peek). The initial and final states of the density matrix 


in the |Dead) /|Alive) basis is: 
5 0 
_ Psys(tp) = i Ls (4.1) 
2 


Psys(0) = ( 


"For the cat, the states |Dead) and |Alive) are stable against further ‘environmental’ disturbance: once 
the cat becomes dead /alive, it remains dead/alive. In addition, (Dead|Alive) = 0 since we take A — 0 limit 
for the cat. 


NIE NIE 
NIRS 
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The simplest lab realization of Schrédinger’s cat is a qubit: |Dead) = | +) and |Alive) = 
| |). However, our ‘system’ is a canonical scalar ¢ defined by (2.8), how do we prepare a 
|Dead Cat) and an |Alive Cat) with ¢? First of all, we need to use pointer states from the 
Hilbert space of ¢, namely states that are least perturbed by the ‘environment’ and therefore 
resemble classical states. There are various definitions for pointer states; and for a generic 
‘environment’, locating the pointer states is a non-trivial problem [1, 38]. Intuitively, the 
pointer states in the present case should be Gaussian wave-packets in @ because they are 
closest to the classical point since AdApg = 2 — the minimal value for a quantum state. 

Therefore we describe the two states |Dead) and |Alive) by two Gaussian wave-packets 


wavefunctions of width o that are centered at ¢ = +p: 
1 _ (¢=49)? f 1 _ ($+)? 
|Dead) = Yı (¢) = a |Alive) = W2(¢) = e 2 . (4.2) 


VN VNR 
The initial state and the density matrix of the cat is then 

1 
v2 


Note that in this representation (Dead|Alive) = 0 only when ġo >> ø, therefore we define 
2 


|Cat) = ¥(4) = =; (%18) + Va(6)) , fsys(0) = [UPEA]. (43) 


$ 
the normalization factor by Trfsysy = 1 and get N = v2ro(1 + eo), The elements of 
Psys(0) are shown in the left panel of Fig. 3. The two ‘diagonal’ peaks represent the two 


wave-packets centered at ¢ = +¢@o (i.e. the high probability of the scalar ¢ being near +¢0) 


o4= o- = ġo => £, A) = (+0,0) (4.4) 
which corresponds to the |Dead)(Dead| and |Alive)(Alive| elements. The two ‘off-diagonal’ 
peaks signify the interference of these two wave-packets: 


$+ = —o- = ġo = (£, A) = (0, £240) (4.5) 


which correspond to the |Dead)(Alive] and |Alive) (Dead| elements. 


4.2 Wigner function: quantum distribution function of phase-space 


One effective way to study quantum decoherence is to go to psys’s phase-space description 
furnished by the Wigner function. So far we have been working with the density matrix 
Psys in the position space: the element p(%, A,t) characterizes the overlap between two 
points (in position space) that are centered in © = Htt- and separated by A = $4 — ¢_. 
Intuitively, when p(%,A 4 0,t) 4 0, the nonzero A signifies the uncertainty in position 
space, therefore should be conjugate to the momentum of the (average) position ©. This 
intuition is borne out by the Wigner function which is defined as the Fourier transform 
(w.r.t. the fast variable A) of psys(X, A, t): 


W(S.p.t) = > I dA eA p, (5, A, t) . (4.6) 


The Wigner function W(%, p, t) is the quantum version of the distribution function f(%, p, t) 
in the phase space {, p} and reduces to f(%, p, t) in the classical limit. However, unlike its 
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classical counterpart f(%,p,t), W(%,p,t) is not always positive-definite — a consequence 
of the quantum uncertainty. Thus, in terms of the Wigner function, one signature of 
the quantum-to-classical transition is the disappearance of its negative part. After this 
process is complete, the Wigner function becomes the probability distribution function of 
the classical states. 

Let us now compute the Wigner function for the cat made of our scalar ¢ — the 
superposition of two Gaussian wave-packets defined in (4.3). In the right panel of Fig. 3 
we plot the Wigner function of the cat at the initial time. This should be compared with 
Psys(0) in the left panel. The two diagonal peaks of fsys(0) defined in (4.4) represent the 
two wave-packets, therefore in the Wigner function they correspond to the two ridges that 


are centered in X = +¢g and extend along p. The two off-diagonal peaks of Psys(0) defined 
in (4.5) correspond to the interference of the two wave-packets, therefore in the Wigner 
function they give rises to the interference pattern along the © = 0 line (note its negative 
parts). 


The cat starts to interact with the ‘environment’ and decohere after t > 0. To see this, 
we first compute fsys(t) starting with Psys(0) given in (4.3), using the convolution (2.14) and 
the propagating function (3.9). Then we Fourier transform to obtain its Wigner function. 
The result has three terms: 


W(%,p, t) = Wı + Wo + We (4.7) 
where 
é 22 l 2 
W(X, p,t) = a wel 7+ 4000) A (a> t) gt a Saza 
AN |h(t)| V rC 
oa | LA E i S2] 57 2 
maag = a lata) a) a] A 


ANA(|VaC |” 


Figure 3. The initial state of a superposition of two Gaussian wave-packets (defined in Eq.(4.3)) 
with ø = 0.1 and ġo = 0.5. Left: Density matrix p(X, A, t = 0); Right: Wigner function W (X, p, t = 
0). 
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correspond to the two individual Gaussian wave-packets and 


2 2 2 
e L 1\- 1 (- =A, - 4, 52-28 Sp- Ap? 
W.(E, p,t) = Auw(t ee ZA ag ii alg 4 (49 
oh +P, ) Aw(t) cos (+E) T Ap] e = ( ) 
with z r 
L $0 1 FA #0 
Ault) a acot ayo? (4.10) 


= e 
AN |h) |V rC 
gives the interference between the two wave-packets. The parameters in (4.8) and (4.9) are 

defined as follows: 


2 . 
C= = fo(0)? + is + aoolt) j (4.11a) 
1. 

L% = A, (4.11) 
L3? = Ze (PAOA) — fa(0)(a10 + a01) — 4Cfi(6)) (4.11c) 
; 4, , 2 ¢ ; 

I = (orc — jolt)? fo(oy? — (ELE q TADAO + a) 
tian , (4.11d) 
and 
penip a (4.116) 
1 = 2C 0 2 0 292 aoo , . 
it = Sf) fr(0) (4.118) 
i} = ~~ 5 (a10 + aor = o? folt) fo(0)) - (4.11h) 


Let us focus on the We term which characterizes the interference between the two wave- 
packets. It has three factors: an enveloping amplitude Aw (t) that deceases over time to zero, 
a cosine factor that produces the oscillation both in X and in p and therefore is responsible 
for the non-positive-definiteness of the Wigner function, and finally an exponential factor 
that encode the shapes and positions of the two wave-packets. To study the decoherence 
process, we only need to watch how the amplitude A,,(t) evolves with time. 


4.3 Purity/Rényi entropy 


One way to characterize the purity of a quantum state is to examine its entanglement 
entropy. The entanglement entropy is zero for pure states and reaches its maximum for 
the completely mixed state. As quantum decoherence is a process from a pure state to a 
mixed state, we may use the entanglement entropy (or the Rényi entropy) to characterize 
it. Then, one may wonder how the entanglement entropy /Rényi entropy behaves when the 
decoherence happens. 
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One hint to answer this question is the analogue with the thermalization, for which 
the thermal entropy saturates when reaching the thermal equilibrium. In a thermal en- 
vironment, quantum decoherence could be closely related to the thermalization, and part 
of the entanglement entropy could be related to the thermal one. Thus, the entanglement 
entropy might saturate when the pure state decoheres in a similar way to the saturation of 
the thermal entropy when reaching the thermal equilibrium. 

There are many possible quantities qualified to be the entanglement entropy, which 
needs to be positive definite and concave. They can all be derived from the reduced density 
matrix. The most commonly used one is the von Neumann entropy: Syn = —Tr sys log Psys- 
The other ones are called the Rényi entropy of order a: 


Sa = i 7108 Tr Êsys (4.12) 
where qa is a positive real number. It can be shown that Syy = limg +1 Sa- 

For simplicity, we will only calculate the Rényi entropy of second order, which is related 
to P := Tepe 3 named as “purity" in the community of quantum information sciences. 
To compare with the results from the phase-space view point discussed in the previous 
subsection, here we also consider the two Gaussian wave-packets for the initial state, i.e., 
(4.3). Using the explicit form (3.9) of psys(t), after lengthy calculations we arrive at 


L2 L2 L? ga, L? g2a 
AT’? La Ly 280-400?) 4 A + F0%00 + 20200 
P = ( ) aie A a ae e 23 + e22 ise 22 Co? a e23 Co? 
LA LA LA LA 49¢2(1+ag9—4C) 
Alt? o*LaLe ,3{ a ai ai ar sd i B (4.13) 
T rN en ee ee E . 
CN? 4 \dll4 — alg 


where the parameters used here are already defined in (4.11). 


4.4 Relation to local quantum quench 


The entanglement entropy has been used to characterize the quantum quench process [39, 
40]. The setup for the quantum quench is to bring the total system to a highly excited 
state and let it evolve. There are different ways to create such a setup. The simplest one is 
to tune a parameter of the Hamiltonian homogeneously such that the original ground state 
turns into non-eigenstate of the new Hamiltonian suddenly. This is the so-called global 
quantum quench. 


“environment” : “environment 


Figure 4. Our setup for quantum decoherence can be also thought as the one for local quantum 
quench in (1 + 1) dimensions. 
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The other is the local quantum quench: Excite the total system locally by either tuning 
a parameter of the Hamiltonian inhomogeneously or simply creating local excitation such 
as quasi-particles or qubits. This setup is the same as ours for quantum decoherence if 
we treat the local excitations (unentangled with the ‘environment’ when created) as the 
‘system’ and the rest of the unperturbed region as the ‘environment’, see Fig. 4. More 
specifically, the initial state of two Gaussian wave-packets given in (4.3) can be thought 
as the local excitations of some effective linear size L which should be determined by the 
profile of the wave-packets, i.e., by o and o. We will try to determine L later. 

One way to characterize the quantum quench is to study the time evolution of the 
entanglement entropy for a region of size L enclosing the initial local excitations or the 
‘system’. By the causal motions of the entangled quasi-particle one expects a cross-over to a 
saturated or mildly growing regime for the growth of the entanglement entropy. The detailed 
time evolution behaviors of the (holographic) entanglement entropy for the local quantum 
quench have been studied extensively [37, 39-42]: the initial growing rate is quadratic in t, 
and then turns to the so-called linear t “Tsunami" regime [40] before dropping suddenly to 
the mildly growing regime. For the case of quantum quench in (1 + 1)-dimensional CFT, 
the entanglement entropy for the a region of size L [37, 40], 


as tx0; 
Sit) ~ Șt, t<L; (4.14) 
const. or Int, t>L. 


For t > L, Sz; either saturates (for global quantum quench) or mildly grows as Int (for local 
quantum quench). Thus, tg ~ L here can be regarded as the relaxation time scale (for the 
probe to relax to its classical state). 

As we can calculate the reduced density matrix for the ‘system’, from it we can then 
evaluate the entanglement entropy for our setup and compare the relaxation time scale tg 
with the decoherence time scale tp. For simplicity, we will consider the Rényi entropy of 
2nd order, i.e., Sg = —logP see (4.12) to extract tg for our decoherence setup. If these 
two time scales are comparable in order of magnitude, then the quantum decoherence is 
closely related to the quantum quench. Though this could be naively expected, it is still 
worthwhile to check this connection directly, especially since there are other environmental 
influences such as temperature which are usually turned off in the consideration of quantum 
quench. We will find that indeed these two time scales are of the same order. Moreover, 


our approach yields similar behaviors to (4.14). 


5 Holographic Quantum Decoherence: Case Studies 


To recapitulate, our ‘system’ is a scalar ¢ whose initial state is a superposition of two 
Gaussian wave-packets (simulating Schrédinger’s cat); the ‘environment’ is a thermal CFT 
with a gravity dual. The decoherence process is captured by the time-evolution of the 
enveloping function A,,(t) of the Wigner function’s interference term and the second order 
Rényi entropy S2(t). The effect from the environment is encoded in the influence functional 
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which is given by ‘environment’s Green’s function. On the one hand, the ‘environment’ is 
non-trivial (in fact strongly interacting), on the other hand, having a gravity dual allows 
its Green’s function to be computed holographically. Once we obtain the Green’s function, 
we can compute A,,(t) via (4.10) and the second order Rényi entropy S2(t) via (4.13), and 
then use them to quantify the decoherence process. 

We will use two examples to illustrate the above. Since the problem is technically 
rather demanding, we present here only simple cases for which analytical computation can 
be carried out as much as possible (in fact up to the point of computing h(7)) and the rest 
is done numerically. Then, to study generic holographic ‘environment’, we also composed 
a python program [43] which can compute numerically h(7) and then A,,(t) and S2(t) for 


any holographic ‘environment’ once its Green’s functions are given. 


5.1 Two Cases 


The two cases we choose are: (1) a static particle coupled to the ‘environment’ of a (1+1)- 
dimensional holographic CFT at finite temperature; (2) a scalar zero mode coupled to a 
(3+1)-dimensional holographic CFT at zero temperature. Their retarded Green’s functions 
were computed holographically and have closed forms, which allow the function h(t) to be 
obtained analytically as well. 


5.1.1 Straight string in BTZ black hole 


Let us first look at a static particle in a (1+1)-dimensional CFT at finite temperature. The 
finite temperature CFT is holographically dual to the ((2+1)-dimensional) BTZ black hole, 
and the point particle is dual to a straight string connecting the black hole horizon and 
the AdS3 boundary.'? Now we summarize the procedure to compute the retarded Green’s 
function holographically. 

The metric of the BTZ black hole is 
(2 


2 p2 
r rH 


pa p2 
ds? = p FT dt? 4 z dX? + dr? (5.1) 
with £ the AdS radius and ry the radius of the horizon. Its Hawking temperature is 


Irl 


(5.2) 


Let us consider a static string configuration connecting the black hole horizon and the AdS3 
boundary. In the static gauge this configuration is given by X(t,r) = 0. Now consider the 
perturbation of X along this static string, the equation of motion of this perturbation 
5X (t,r) can be obtained from the Nambu-Goto action of the probe string [35]: 


2 2 
-pp EH at Oy (1? (r? — ri7)Oy) | 5X (tr) = 0. (5.3) 


12This setup was also used in [35, 36] to study the holographic quantum Brownian motion (which has 
been briefly reviewed in Sec. 3.1). The horizon-end of the string receives random kicks from the Hawking 
radiation near the horizon, and these random kicks then propagate along the string and cause the Brownian 
motion of the boundary-end of the string. 
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The solution with the in-falling boundary condition at the black hole horizon is 


iw? 
=i f r—iwl? (r—ry\ 2H 
5X (t,r) = e“ X,,(r), with Xl) x m (2 = =) . (5.4) 
The conjugate momentum of Xw (r) is 
Po o e 
IL, (r) = -Ts l n) d,X.,(r) (5.5) 


where T; = 2 is the string tension. Then using the prescription given by (2.32), we obtain 
the retarded Green’s function: 


. I(r) wre + iwrr 
G =- 1 = c g 
RCs re) re wt) * riwl 
r2 5 9 r2 
= Tre ( 14 2) w° — iT, ( wer ae u) (5.6) 
Cc 


with re (œ~ oo) the UV-cutoff. 

The real part of (5.6) is UV-divergent but can be considered as the renormalization of 
the probe’s mass, and only the imaginary part of (5.6) drives the quantum decoherence. 
However, its form has unreasonable high frequency behaviors and we need to regularize by 
introducing some window function, which we choose to be the commonly used Lorentzian 
function (with width T,,).1° With the real part dropped and after regularization, the 
retarded Green’s function (5.6) becomes: 

N2 (, 3 22 r3, 

Gr(w) = -i N4 (w? +4r’°T’w) Te +2 (5.7) 
where Ng = VT; = ¢/¢, is the number of degrees of freedom of the dual CFT. Then the 
symmetric Green’s function Gsym can be obtained from (5.7) via (2.30). 

With the closed form of Gr given by (5.7) we can solve for h(7) via the inverse Laplace 
transform (3.13). First, the Laplace transform of G(r) is 


z 1 
Gr(s) = aN202, (40°T? -T? . 
R(S) TIN st w ( T w) Pu +s 


We then compute h(7) using (3.13). (Recall that once h(T) is known, all other functions 


(5.8) 


can be immediately expressed in terms of simple functions of h(7) and its derivatives.) Let 


us consider the case given for which h(s) has two complex poles at s = sı +is2 and one real 
pole at s = s3, where s;’s are real functions of T and I. However, their detailed expressions 
are not very illuminating and will be omitted. Then the inverse Laplace transform of h(s) 
gives A(T): 


so(s3 +Tw)(e%7 — e817 cos s27) + (8? + 8% — s3l'w + 51(T'w — 83))e%!7 sin soT 


s2(8? + s2 — 25183 + 52) 


h(t) = 


(5.9) 
Knowing h(r), we can then first compute f;’s and gis by (3.14), and then use them to 
calculate the Wigner function and the second order Rényi entropy. 


13These unreasonable UV behaviors are termed ohmic (~ w) or super-ohmic (~ w? with 5 > 1) in the 
literature, see e.g.[8] which also used the window function to regularize the UV-behavior. 
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5.1.2 Scalar probe in AdS; spacetime 


Most of the earlier studies of quantum decoherence assumed that the ‘environment’ is a 
thermal reservoir. This might cause the misconception that the thermal fluctuation of the 
‘environment’ is necessary for the quantum decoherence of the ‘system’. However, as we 
are about to see, decoherence can happen even in zero temperature. 

To show this, we choose the ‘environment’ to be a (3 + 1)-dimensional CFT at zero 
temperature and the ‘system’ a scalar operator. The total setup is holographically dual 
to a scalar in AdS; spacetime. Following the procedure described in Sec. 2.3, one can 
compute holographically the retarded Green’s function GR of an operator O with conformal 
dimension Ag. This was done in [28], and we summarize the results here: 


N20 (3-Ag)e2(40~4) (q2)40-? [ cos TA — i sgn(w) sin TAO ] 2 < Ao ¢ N 


Ga(u) = ¢ EGO 
a) Nz So) (q2)40-? [ Ing? — i r sgn(w) ] 2<AoeN 
(5.10) 


812 (Ao—3)!22740-5 
where N2. is the number of degrees of freedom of the dual CFT, and € ~ 0 is the UV cutoff. 
For simplicity, we only consider the zero-mode (with q = Vw? — k2 = w) of ¢, namely 
we assume the ‘system’ is homogenous in space. The problem is then effectively (1 + 1)- 


dimensional, which can be compared directly with the straight string case. In contrast to 
the static string case, the non-zero and non-renormalizable real part of the retarded Green’s 
function is absent here. 

In order to solve A(T) analytically by the Laplace transform (3.13), we only consider 
the integer and half-integer Ao cases. Moreover, as Gr contains the ‘super-ohmic’ factor 
(gre? we will again regularize by introducing the Lorentzian window function of width 
Tw. The Laplace transform of G(r) is the following: 


N20 (3-Ao)FR40—Y pant+2 1 1 
Gr(s) = SnD (Ap -2) 250-5 lw sFr 2<AoeN+o 
B N2,€240-4) Aop2Ao-3 2 1 1 1 
s2403 22505 1)Poly mht, ar, tar, ar) *S40eEN 


(5.11) 
Then h(t) can be computed by inverse Laplace transform (3.13). For large Ns- and Dw, 


h(s) has two conjugated poles at s = sı + is2 and two real poles at s = s3, 54, therefore 
h(t) takes the following form: 


[s2(s3 — To) (52 + (s1 — 84)” )e® + s2(83 + (s1 — 83)°)(T%, — sz)e*” 
+ (st — s1(s3 + s4)(s? — s2 + T2) + (s2 + T2,)(s2 — 5384) + s?(28% — T2, + 8384)) (s3 — s4)e* sin sot 
(s{(s3 + s4) + (85 +T2,)(s3 + 4) — 281(T'2, + 5384)) 82(s3 — 54)e®™ cos sot | . (5.12) 


Again, we omit the detailed expressions for s;’s. In the numerical calculation, we set 
Nee = 10 and the UV-cutoff € ~ 1/Tw. 
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Figure 5. The evolution of the Wigner function W(X, p,t) for the two Gaussian wave-packets 
made of the scalar ¢ with Ao = 3. The evolution goes clockwise as t = 10, t = 15, t = 20 and 
t = 30 units respectively. 


5.2 Results and remarks 


Based on the the expressions of h(t) in the previous subsection, we can now calculate the 
enveloping function A,,(t) via (4.10) and the second order Rényi entropy S2(t) = — log P via 
(4.13), and use them to characterize the decoherence of the superposition of two Gaussian 
wave-packets. 

The task requires long numerical calculations. Before we present the detailed results 
below, let us first show (in Fig. 5) four snap shots in the time-evolution of the Wigner 
function. Note the decay of its negative part (which characterizes the quantumness of the 
Wigner function). In the scenario of the environment-induced decoherence, the coherence 
decay exponentially. Therefore, we first need to define precisely what we mean by deco- 
herence. Since the decay is exponential, only the order of magnitude of tp is relevant. We 
will declare that the ‘system’ has almost decohered once the enveloping function A, (t) 
decreases to the 1% of its initial value and define the decoherence time tp by: 14 


Aw (tp) 
Aw(0) 


decoherence time tp : = 0.01. (5.13) 


Note that the time-dependence of the Rényi entropy S2 exhibits a crossover around t % tp. 

From our computation we extract two important results. First, for the decoherence 
time tp, we obtain its dependence on the temperature T of the ‘environment’ and on the 
conformal dimension Ao of the field © coupled to the ‘system’. Second, we obtain the 
growth rate of S2(t), using which we can define different stages of decoherence. We then 
show that these stages are nicely matched with those in the holographic quantum quench 
discussed recently in [40-42]. 


“The decoherence time tp is roughly of the same order as the half-life scale of Aw(t) since Aw(t) is 
exponentially decaying. Here defining the ‘almost decoherence’ by (5.13) is simply for convenience. 
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Figure 6. Time evolutions of the enveloping function A, and the second order Rényi entropy S2 
for a string probe in BTZ background. The parameters chosen for (4.10), (4.13) and (4.12) are 
(60,0, Net, Q, Pw) = (05, 0.1, 1,1, 1,10). Left: A,, (thick lines) and S> (thin lines) for T=50 (circle), 
30 (triangle) and 0.1 (square). Note that we have rescaled S2 by 1/2. Right: The long-term time 
behavior of S2 for T =20 and 1, which will be used for the fitting in Fig. 8. 
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Figure 7. Left: Decoherence time tp v.s. temperature T. Right: Zoom-in of high -T regime 
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Note that in all the numerical plots presented below, we set the width of the window 
function T'ẹ = 10, the coupling g = 1 and the constants Ng = 1 and N.e = 10.15 The 


width I’, also serves as a reference inverse time scale. 


5.2.1 String probe 


For the case in which the ‘system’ plus ‘environment’ are the holographic dual of a string 
probe in a BTZ background, we present in Fig. 6 the time-evolution of the A,,(t) and S(t) 
at different temperatures. We can see that Aw almost decays away at t = tp, and around 
this time scale S2 goes through a crossover. As discussed, this implies that the quantum 
quench, which is characterized by this crossover behavior, is closely related to the quantum 
decoherence. 


Next we study the temperature (T) dependence of the decoherence time tp. The results 
are shown in Fig. 7. First, let us look at the left panel which covers the entire 0 < T < 200 


Note that once we have taken the ‘weak coupling limit’ by keeping only the quadratic terms in nF 
(and reaching (2.21)), we are free to absorb the coupling constant g into the scalar fields, i.e. we can simply 
set g = 1 for convenience. 
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region. We see tp decreases as T increases, namely, the hotter the ‘environment’ the faster 
the ‘system’ decoheres — consistent with the intuitive picture of decoherence. 

Next, we look at the high-T and low-T regimes separately. A conformal field theory 
has no intrinsic scale which we can use to define high-T or low-T regime. However, when 
we regularize the Green’s function with a window function, we introduce a length scale into 
the system in term of its width Tw; therefore we can define high-T or low-T w.r.t. Tw. 

In the low-T (T < Tw) regime, we can see the decoherence time tp is almost independent 
of T. The reason is the following. The window function of width I, introduces a ‘soft mass 
gap’ into the CFT. When T < Pe the excitation is suppressed by the mass gap and the 
decoherence process resembles the one at zero-temperature therefore is insensitive to the 
actual temperature. 


For the high-T regime (see right panel of Fig. 7) we discover a nice scaling behavior 
1 


by fitting the data points. This scaling behavior is the same as the one extracted directly 
from the master equation (3.16) for the simple toy model in [1]. While not surprising, this 
is the first time this scaling behavior is obtained when the ‘environment’ is a non-trivial 


CFT. 


Finally, let us look at the time-evolution of the second order Rényi entropy S2 of the 
‘system’. In Fig. 6 we can already see that S2 has different scaling behaviors at different 
regions of t. This suggests that the decoherence process happens in stages. With the fitting 
tool ‘formulize’ [44], we fit the data and divide the decoherence process into the following 
four different stages according to the scaling behavior of S2(t): 


Co ? + O(t’) , t< 2 (initial) ; 
t 
S(t) ~ Coc. t+ O(1), ists (intermediate) ; (5.15) 
Cac. nt+O(1), tp StS 10tp (late) ; 
slower than C Int , 10tp <t (final) . 


where Co, Co.c., Cae. Ci(< Cac.) are all positive and increase monotonically with T. The 
data (together with their fittings) of the first three stages are plotted in Fig. 8. We also 
compare the high-T and low-T cases and find that the scaling behavior of S(t) is insensitive 
to the temperature. The intermediate stage is about one order of magnitude longer than 
the initial stage, which agrees with the corresponding result for the quantum quench in 
[40]. After the crossover until the end of the late stage (at ~ 10tp), S2 grows as Int, in 
agreement with the results in [37, 41]. During the final stage, we see the log t-growth of the 
previous stage slows down, although our present limit with numerical computation (due to 
accumulated numerical errors that affect particularly the large-t region) prevents us from 
determining the long-term trend of S2(t).! Matching the first three stages of (5.15) with 


16Tn particular, S2 acquires a negative linear-t term to offset the positive Int term. However, instead of 
the eventual saturation that we (albeit maybe naively) expect, we see a Cot — C3 Int growth at the very end 
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Figure 8. Scaling behaviors of S at different stages of the growth shown in the right panel of 
Fig. 6. The three rows from up to down corresponds to the three stages defined in the main text. 
Initial: The solid (orange) lines corresponds to $2 ~ t? + O(t*) and fit better than the solid (blue) 
lines that have only t? term. Intermediate: S2 ~ t + O(1). Late: S2 ~ Int + O(1). 


(4.14) we conclude that the relaxation following a local quantum quench occurs when the 


local excitations decohere. 


5.2.2 Scalar probe 


Similar to the string probe case, we first present in Fig. 9 the evolution of A,,(t) and S(t) 
for operator O with different conformal dimensions. The qualitative behavior is the same 
as the string probe case. Again, we find that the quench and decoherence time scales are 
of the same order. In fact, since now the ‘environment’ is at zero temperature, it is closer 


(using 20 additional data points extending the right figure of Fig. 6). The two terms Cot and C3 lnt are 
roughly equal in magnitude and have opposite signs therefore we cannot predict the long-term trend based 
on this very end; and although this could be a mirage caused by the large-t numerical errors, we mention 
this as an issue that deserves future study. 
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Figure 9. Time evolutions of the enveloping function A, and the second order Rényi entropy S2 
for a scalar probe (zero-mode only) in AdS; space-time. The parameters chosen for (4.10), (4.13) 
and (4.12) are (¢9,0, Nsc, Q, Tw) = (05,0.1,10,1,1,10). Left: A,, (thick lines) and S (thin lines) 
for Ao =4.5 (triangle), 4 (circle) and 3.5 (square). Note we have rescaled S2 by 1/2. Right: The 
long-term behaviors of S2 for Ao =3.5 and 3, which will be used for the fitting in Fig. 11. Note 
that we have rescaled the time for Ao = 3 case by 1/10 in order to fit it in the same plot with 
Ao = 3.5 case. 


Figure 10. Intp v.s. Ao for 2.5 < Ao < 5. The blue line is the fitting function: In tp = 
2.26 — 3.3440. 


to the quantum quench setup. This provides more supporting evidence for the connection 
between quantum decoherence and quantum quench. 
Next, we fit the scaling behavior of tp w.r.t. Ao and find (as shown in Fig. 10): 


Intp ~ — Cado (5.16) 


with Ca > 0. (This is to be compared with the scaling tp ~ 1/T in the previous case 
(shown in Fig. 7).) This means that the pure state of the superposition of two Gaussian 
wave-packets has a lifetime of order A Namely, if we try to realize a qubit embedded in 
strongly-coupled environments of the type considered here, we should choose those with as 
small a Ao as possible, in order for the ‘system’ to be more robust (against decoherence). 

Finally, using the fitting tool ‘formulize’ we determine the scaling behaviors of the 
second order Rényi entropy S(t) at different stages for two different Ag’s and the results 


are shown in Fig. 11. Although the total system is different from the previous case, we 
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Figure 11. Scaling behaviors of S at different stages of the growth shown in the right panel of 
Fig. 9. The three rows from up to down correspond to the three stages defined in the main text. 
Initial: The solid (orange) lines correspond to S> ~ t? + O(t*) and fit better than the solid (blue) 


lines that have only t? term. Intermediate: S2 ~ t+ O(1). Late: S2 ~ Int + O(1). 


found the scaling behavior to be basically the same as (5.15). 
C’s in (5.15) increase with Ao.) This suggests that the scaling behaviors of S and the 
definition of the different stages given in (5.15) might be universal for different probes and 
environments.!’ It is also very reassuring to see all our results agree with the studies based 


(And all the constants 


on the (holographic) entanglement entropy for the local quantum quench [37, 39-42]. It 


would be interesting to explore this connection further. 


'’The lnt growing behavior in (4.14) or (5.15) could be a feature of 2D CFT. Although we consider the 
scalar probe in AdS; at very beginning, we only focus on the decoherence behavior of its zero-mode, which 
can be thought as an effective (1 + 1)-dimensional problem. Thus, we obtain the scaling behaviors (5.15) 
as expected. 
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5.3 Initial state profile dependence of decoherence time 


The decoherence time tp depends on how the ‘system’ (the scalar ¢) interacts with the 
‘environment’ and on the properties of @ such as its mass or conformal dimension. The 
other relevant attribute of the ‘system’ is the profile of its initial state, e.g. in our case 
of two Gaussian wave-packets the relevant data is the width (øo) of each packet and the 
‘distance’ (269) between the two. 

Further, the similarity between (5.15) and (4.14) implies a close relation between the 
quantum decoherence and the local quench once we identify the decoherence time tp in 
(5.15) with the effective length L in (4.14). Then once we determine tp’s dependence on 
the initial state profile, it immediately gives the relation between the initial state profile 
and the effective length. This dependence is nontrivial due to the complexity of A,,(t) given 
in (4.10), and needs to be extracted by numerical methods. 

In Fig. 12 we show the decoherence time tp as a function of the profile ({o, ¢o}) for 
the scalar probe with Ag = 3. Fitting the data by constant ø or ġo slices, we find 


c(a) 


po + c2(0) on G D ~ ¢3(¢0) + ca($o)o or fixed ¢o. (5.17) 


tp 
This shows that the effective length (+ tp) is not proportional to ġo (the distance between 
two packets) as naively expected for uniform probe. Instead, the relevant quantum infor- 
mation is only encoded in the properties of each individual Gaussian wave-packet, and the 
decoherence time tp or the effective length L increase with ø. Thus, for fixed o a larger 
oo means a bigger zero-information region therefore lower concentration of nontrivial in- 
formation, which accelerate the decoherence or the quench process. Therefore, to suppress 
decoherence, we should prepare initial states with small ¢o. 


0,0006 


Figure 12. Dependence of decoherence time tp on the profile ({a, do}) of initial state of the two 
Gaussian wave-packets, for scalar probe with Ag = 3. All other parameters (except for o and ¢o) 
are the same as used in Fig. 9. 
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5.4 Numerical studies for discrete time systems 


All the above results are based on exact solutions of h(7) without any (numerical) approxi- 
mation. This is possible because the retarded Green’s function GR is simple enough to allow 
h(r) to be computed analytically via Laplace transform. In general, Gp is more complicated 
and solving (3.13) to obtain h(T) requires numerical methods. In order to generalize our 
proposed scheme in this paper and examine the environment-induced decoherence scenario 
for more general situation, we compose a python code [43] to solve (3.13) and then evaluate 
Aw(t) and S(t) for a given retarded Green’s function. 


10 


Wigner_Renyi2 
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t 


Figure 13. Evolution of Aw (decreasing) and S2 (increasing) for the two Gaussian wave packet of 
scalar probe for Ao = 10/3. The discretized time unit is 1075 and Fẹ = 4000. 


Here we show a typical example computed by our python code [43]: the retarded 
Green’s function is given by (5.10) with Ag = #. Note that for this case (3.13) can only 
be solved numerically. The results for Ay(t) and S9(t) are shown in Fig. 13. It again 
shows the decoherence behavior as expected, however, the detailed behavior depends on 
the discretized time unit and I,,,.!° 

The python code still needs to be improved, in particular, its sensitivity to the choice 
of the (discretized) time unit and of the window function needs to be further reduced. Once 
this is achieved, we could use it to search for the robust systems against generic quantum 
decoherence — a necessary ingredient in designing real quantum computers. 


6 Conclusion 


In this paper, based on the earlier observation |23, 26] we elucidate the connection between 
Feynman- Vernon and Schwinger-Keldysh formalisms: The influence functional of Feynman- 
Vernon is the generating function of the Schwinger-Keldysh correlation functions. We 


18Despite the efficiency of this numerical code, we still need to discretize the time (with the main bottle- 
neck being the integral-differential equation (3.12)), which leads to a discretized retarded Green’s function. 
Therefore, the final results depend on the unit scale of the discretized time and would be different from the 
continuum time limit. 
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then apply this connection to study quantum decoherence when the ‘environment’ is a 
holographic CFT, since the Schwinger-Keldysh retarded Green’s function can be obtained 
holographically. We consider two cases for which the influence functional can be computed 
exactly, and derive the time evolution of the ‘system”’s reduced density matrix, from which 
we then evaluate its Wigner function and Rényi entropy. We use the disappearance of the 
negative part of Wigner function to characterize the decoherence process and define the 
decoherence time. We then extract the scaling behaviors of the decoherence time scale. 
Although we focus on these two particular examples, we expect the results to be rather 
universal for strongly-coupled CFT ‘environment’. The second example also proves that 
the quantum decoherence occurs even at zero temperature. 

The ability to determine the dynamics of the reduced density matrix allows us to di- 
rectly compute the Rényi entropy which was used to characterize the quantum quench. 
As explained earlier, the setup for the quantum decoherence considered in this paper is 
very similar to that of the local quantum quench. We demonstrate the close relation be- 
tween these two quantum processes by comparing their characteristic time scales and the 
scaling behaviors (with time) of their Rényi entropies. Our results indicate that the lo- 
cal quantum quench occurs when the local excitations decohere. Moreover, we show that, 
for the quantum decoherence process, the long term scaling behavior of the Rényi entropy 
deviates from the logarithmic growth exhibited by the local quantum quench process. It 
would be interesting to further explore the relation between these two dynamical quantum 
phenomena. 

Now a few words on the relevance of the present work to the design of realistic quantum 
computers. Although there are many different lab realizations of qubits, such as iron traps, 
liquid NMR, optical lattices, to name just a few, in order to have large number of qubits in 
the ‘system’, solid-state based designs, such as semi-conductors and superconductors, are 
more promising candidates. Among solid-state structures in the lab, many are strongly- 
coupled and have (at least approximated) holographic duals. The method introduced in 
this paper can be applied to study the quantum decoherence in any strongly-coupled ‘en- 
vironment’ that has a holographic dual; therefore it is useful in the systematic study of the 
general behaviors of quantum decoherence and could hopefully provide some guidance in 
selecting the type of environments that are more coherence-friendly. 

The most interesting open problem is how to incorporate the probe holographically. 
Answering this question would also solve the following technical problem. In this paper, 
we needed to implement numerical computations to determine the time evolution of the 
reduced density matrix, which require solving integral-differential equations and usually 
cause numerical artifacts. To bypass this technical difficulty, one could try to find a method 
to directly compute the propagating function (2.14) holographically, which is possible once 
we know how to describe the probe ¢ holographically. We hope to explore this possibility 
in the near future. 
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A Evaluating the propagating function: a review 


In this appendix, we evaluate the propagating function J by regarding it as a partition 
function (3.3) with action (3.4). This has been done before, e.g. in |7, 8]. Now we review 
the derivation to make the paper self-contained and to fix the notation. 

First, we go to the momentum space, in which the effective action (3.4) becomes 


t 
Sett -5 | dr [;A;- (9? — |kP?)DzAz 
i 


t t >, ; 
-¢ >| ar f dr farina er = TERT) = sAR(TIGQ(T = TART!) 
i 


(A.1) 
As explained in the main text, now we will only consider the k = 0 mode to focus our 
attention on the time-dependent behavior. Then evaluating (3.3) only takes two steps: first 
we expand the action (A.1) around the classical saddle points, then we integrate out the 
fluctuations to obtain the overall normalization factor. Expanding © and A as follows: 


y= +m, A(r) =A+q@ (A.2) 


where X and A are the classical solutions that solve the EOMs derived from the real part 
of (A.1): 


(7) +9? ET) +g? a dr'Gr(r — TET) =0, (A.3) 
A(r) +0? A(r) +9? i dr'Gr(7’ — T)A(7’) = 0 (A.4) 


with the boundary conditions (in time) 
p p (A.5) 


The time variable r € [0,t]. Note that for the ease of reading, in the derivation of this 
Appendix we use {Xo,1, Ao,1} but will switch to the tilde and bar notation in the main text. 

Now let us solve the boundary value problem (A.3)-(A.5). First, note that (A.3) and 
(A.4) are related by a change of variable r — t — r, therefore we only need to solve one of 
them. The matching of the boundary data can be done as follows. First we separate the 
boundary data out of {2, A} by introducing { fo,1, 90,1}: 


X(T) = fo(r)Bo + fi(7)¥1 and A(T) = go(r)Ao + gi(7) Ai (A.6) 

then { fo,1} are two solutions of (A.3) and {90,1} of (A.4), and with the boundary conditions 
fo(0)=1, fot) = 90; f0) =0, fi@j=1, (A.7) 

90(0)=0, — go(t)=1; g(0)=1, g(t)=0. (A.8) 


— 36 — 


Combining (A.3), (A.4) with (A.7) and (A.8), we see 
go(t)=fitt-7) and = gilt) = fo(t—T). (A.9) 


Therefore we only need to solve for f;’s. 
fis satisfy the integral equation 


f+ fir +9 f "dr’Galr —1')f(7!) =0 (A.10) 


with initial and final values given by (A.7). However a solution to (A.10) can be uniquely 
determined either by fixing the initial and final values (boundary value problem) or the 
initial value plus its derivative (initial value problem); since a boundary value problem 
usually requires the shooting method therefore is rather difficult, we recast this boundary 
value problem into an initial value problem (by shifting the information at T = t to T = 0) 
following |7]. 

To convert (A.10) into an initial value problem, let’s first apply the Laplace transform 
on it, suggested by the convolution term. The Laplace transform f = L(f) of f satisfies 


i 1 g 
r= 52 + 0? + g?Gr(s) ea s? + 0? + g?Ga(s) a 


The original boundary value problem is now recast into an initial value problem: for given 


(A.11) 


f(0) and f(0), f(T) is uniquely determined as the inverse Laplace transform of (A.11). It 
has two fundamental solutions: 


h(0)=0 A(0)=1 = h(r) = L~ a (A.12) 
H(0)=1 HH(0)=0 = H(r)=£ (sama 


We see H(r) = h(r), since H(s) = sh(s) and h(0) = 0. 

Since the integral equation (A.10) is linear, all its solutions are linear combinations of 
h(t) and h(r). The two solutions to the boundary value problem (A.7) are therefore easily 
determined to be: 


h(t) 1 
2 h = — h A.1 
aD =O n(n) + fiir), sale) = GA) (A.13) 
where h(rT) is the solution to the initial value problem 
h(r) +9 h(r) +e | dr'Ga(t —7')h(7') =0, h(0) =0, h(0) =1. (A.14) 
0 
; TEO , : — 1 
For simple Gp, (A.14) can be easily solved via inverse Laplace transform h(T) = £ (saree ) ; 


Note that due to the boundary condition (A.7), {fo, fi} are actually functions of both 7 
and t, with t entering as the boundary point. 

Substituting (A.2) into (A.1) (with k = 0 mode being considered) and use the equations 
of motion (A.3) and (A.4), one obtain 


Set = Sa[X, A] + Salgı, q2] (A.15) 
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with 
g i t t 
Sale, Al] = S(rIA(r)|s + sv f ar | dr’ A(7)Gsym(T —7')A(1’) (A.16) 
0 0 
and 


Salata] = far [dar ir) - Palet) 


a 


= | ar | ar faaleyGntr -Da = ionlae) 


(A.17) 
Note also that one can obtain the same results by using the complex saddle point. 
One can factor out the boundary data in Sea and obtain 
1 st 
Sa = >_[filt)SeA1 — fi(0)ZiAo] + A Soo aAA; (A.18) 
i=0 i=0 j=0 
with 
t t 
aij(t) = Pf ar | dT'gilT)Gsym(T — F )g3(7’) . (A.19) 
0 0 


On the other hand, to carry out the integration over the fluctuation, it is convenient to 
rewrite Sq into the following form 


) EPS. T ‘ar! Tir TaT 7! 
isa= -3 f dr | a QM OC -7)QC") (A.20) 


where QT (t) = (a(t), qı (t) ) and 


oeie 9? Geym(7 — 1!) i[(02 + PJA- 1!) + GRET — 7’) 
i[(OF + 07)5(7 — 7’) + 9°Gr(T —7')| 0 l 
(A.21) 
In the Appendix C, we show that the normalization factor 
; 1 
At) = | Dq Dq eal”) = =1/2 as A.22 


with h(t) given in (A.12). 


Combining the above results, we obtain the propagating function J 
J(X1, A1, t|£o, Ag, 0] = A(t) ar Diol fi(t)EsA1—fi (0): Ao]—$ Zizo Zj-o aiz (t)A:4; (A.23) 


with a,j(t) given in (A.19) and A(t) in (A.22). 

In summary, the propagating function (A.23) can be obtained by just solving h(t) and 
H(r) for a given Gr using (A.14). Then, the time evolution of the reduced density matrix 
is given by (2.14). 
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B Derivation of the master equation for /sys(t) 


Now we derive the master equation of Psys(t). Start with the propagating function (A.23) 
and recall that the function f;(t) in (A.23) is actually a shorthand for fi(t) = filr = 
t,t) where { fo(7,t), fi(7,t)} are the function { fo, fi} given in (A.13): fi(7,t) = a and 
fo(7,t) = h(r)- mht) and the dot refers to the derivative w.r.t. the variable 7. Therefore 
we have 


h(t) 1 
AOE Aha Ë, TES 
‘ [H(t)|2 ) (B.1) 
fo(t) = H(t) — nt)’ fo(0) = — e= ft E 


Now evaluate the t-derivative of the propagating function (A.23): 


. d L; Al iSeff 
ill) =iA (te 


+i ~ (t) + AZo fo(t) — AoE fi (0) — iaso Zilo) J(t) (B2) 


T i(-3 = Ajai (t = A Ag oW ce z AoMo(t )) J(t) 


where the prime denotes t-derivative, with 


Aw) = ZA = 4 EA =AW-A@P, AQ) = -AWAL0), 
LO a + folt) fo(0) , f4(0) = —fo(t) fo(0) , (B.3) 
A't) = it = —f,(t)A(t) . 


Recall that f;(7,t) satisfy the equation of motion 


d Q2 2 P G _ BA 
Soin) +P Flr) + 9 i saeg oi (B.4a) 


with boundary condition fı(r = 0,t) = folt = t,t) = 0 and Air = t,t) = folT=0,t)=1. 
Using this we compute the double derivatives in (B.3): 


2, t 
A= LAO E e / Pen eT AC eee Cee 


d2 


; (B.4b) 
jolt) = Fah) = -9 f ds Grlt— sols) = bol 


where 0;(t) = g 2 ds Gr(t — 5) fi(s). 
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Now let us derive some useful formulae: 


ap) = [ifi(t)Aa — ifi (0)Ao]J , (B.5a) 
a = BAO + i fo(t)Xo —ayzA1 — SoTa Ag 3 (B.5b) 
1 
x sd = HAH- [AOPA — fit) fot) Aro + fit) fi (0)Z1 Ao + folt) fi(0) Zo Ao] J 
1001 
+[iarrA1(fi(0)Ao — fi(t)A1) + 5 (a10 + ao1)Ao(fi(0)Ao — AHADI 
= [ifi(t) — AOPA — fi) fo(t)A1Do + AOAO Ao + folt) fi(0)ZoAo] J 
-É (aro + ao1)Ao ~ + ayAy = |J. (B.5c) 
With (B.3), (B.4) and (B.5c), the first two lines in (B.2) can be reduced into 
a JANA A1 £1 J+b1 Aad bp hp AiJ— [= z (a10-+001) Ao +a di cold P (B.6) 


Then using 


dal) = 9S f fasts oa(9)Gaym (580218) = 280l) f Gay (ts)oa(s) =0 (BT 


and (B.5a), we can rewrite the last line in (B.2) into 


(Geis Set HO) ots, SP oi 8 5. (B.8) 


7 | z 1 
fO O 2f (0) 02 
Finally, combining (B.6) and (B.8) and using (B.5) to replace all No and Ao, we obtain 
d ð ð . bo o 
<J = Q? +IP AN S J 4 : A J 
Sea yan] r+ ( a) ‘BA, 
å10 + Gor ae 2 filt) bo(a10 + a01) r ð J 
2 fı (0) fiO) (aio +ao1) 2f (0) folt) Oxy 


(boa , aofolt) 1. | dio +401 filt) A27 
l (oe “AAO 277 2 20) ! 
¡210 + %1 o? 

2fi(0) Ot 


(B.9) 


where 6? = bı — bo ne Then it is straightforward to write down the master equation for 
the reduced density matrix psys(X1, A1, t) from (B.9) and (2.14). 


C Normalization of sys 


C.1 Evaluating the normalization A(t) using Trfsys(t) = 1 


We now use the propagating functional given by (A.23) and the condition Trfsys(t) = 1 
to evaluate the normalization factor A(t). Here we assume the initial density matrix is a 
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Gaussian wave packet. 


n : 1 _ $4482 1 _ Agtanz 
(b+ Peys(0)|@_) = =e" z “cr (C.1) 


O/T ~ N 


where the normalization factor ø is fixed by the unitarity condition Trfsys(0) = 1. 


The reduced density matrix is 


($+|Asys(t)|o-) = f dXo dâo J[E1, A1, t|Z0, Ao, 0} (6+ lôsys(0)lġ-) 


Ae 2 Ath (t)A1E1— 2 fo(t)2A? 


(C.2) 


J dage (FAO) AF (F fo(t) fo(0)Ar— gar0lt) aos Aif 0) Ao 


TOTI A(t)e7 tan (t)A2+i fi (t)A1D1— oF fo(t)?A2 


(if: (O)S1 + F fo(t)fol0)Ar — aoli) + aor(#))A1) 
AC’ 


xX exp 


il y2 a YA 1 A2 
141 
Le, 1 L2 1 


= TO? A(t)e ia 
where the characteristic length Ly, Le, La and C are defined in (4.11). 
Now using Trfsys(t)=1, we can derive the normalization factor A(t). 
1=Trpys(t) = f dodd < dlpay(t)|6! > IO- A), (C.4) 
= f manista) < ÎÊsys(t) >, 
‘ 2 
—if(0) 
= f vrane. exp (A0) a ; 
AC 
= Á. 2r — : f 
| f1(0)| 
Using (B.1), we have 
Al) = = (C5) 
mhk, 


C.2 A consistency check: Evaluating the normalization through functional 
determinant 


In this subsection, we compute the normalization factor A(t) by evaluating the functional 
determinant in (A.22) directly. We follow the formulation in [45]. 


Suppose we want to compute the determinant of an operator L 


Lao; + P(t), with t€ [to, ty] (C.6) 


— 4] — 


where P(t) is a real function and the boundary condition for the eigenfunction of L is 
u(to) = u(t) = 0. The prescription takes only two steps. 
First we write the boundary condition as 


u(to) ult) \ _ (0 ' 
m (a) av (| = (4 em 


where M and N are 2 x 2 matrices and in our case 


M= Gy and N = Gi ; (C.7b) 
00 10 


Then we define matrices H(t) and Y(t f) as 
7 yi(t) y2(t) > FiA 
H(t)=|~ Í i Y(t) = H(ts)H (to) (C.7c) 
(’ (t) ya(t) i i 
where y1 (t) and y2(t) are two independent solutions of the homogeneous differential equation 
L{ ™ | =o. 
y2 
The determinant det Ê is certainly divergent, so the only meaningful quantity is the 


renormalized determinant, i.e. the ratio (C.7d). Usually we chose L = 8? as in [46]. Finally 
we have the determinant ratio of È and some other (similar type of) operator L: 


det Ê  det[M + NY(ts)| 


== = ; C.7d 
detL det|M + NY (t+)] ( ) 
Here the matrices M and N are given in (C.7b) for our case. The determinant is 
^ t tr) — yalt t 
det[M + NY(t,)] = yı (to)y2(tz) — y2(to)yı (ts) (C.8) 


~-yr(to)¥2(to) — y2(to)ýı (to) ` 
If we pick y;(t) being the solution whose initial condition is yı(to) = 0, the determinant 
ratio (C.7d) becomes 


dettL 1 y(tr) 


— = : C.9 
det L det L y (to) oe) 


Now the normalization factor obtained by integrating out the quadratic action (A.20) 


- 1 
A(t) = [det(0)] 7" = dal 5G =) 4 Ga OI (C.10) 


with O(t — t’) defined in (A.21). 
It is obvious that the solution y(t) in (C.8) is given by the solution f(t) in (A.13). 
Hence the normalization factor (C.10) becomes 


. =l 
a AI (1O) gri =i 
a IAO co 7 Eo (C.11) 


which agrees with (C.5). 
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C.3 Another consistency check: The behavior of J(t) and fgys(t) as t + 0 


As a second consistency check, we compute the behavior of the propagating function at 


t — 0 and check that it indeed recovers a Dirac delta function. Using (B.1) we have 


: ; 1 : : 1 
dim f1(0, ¢)| = € i Jim | filt, t) = < , 
. 1 . 1 
i ——— l = —— .12 
Jim | fo(0, t) € ’ Jim | folt, t) č , (C ) 
: 1 
oe A(t) 7 Ine f 


and aij(t + 0) + 0. The propagatinf functional (A.23) becomes 
PEA ei(2(A1—Ao)(=1—Xo)) ss eit (O4+-$4)? . e-it(o--4-) 
t=e0 QTE t=30 = \/2rie /2n(—i)e (C.13) 
= 6(b4 — $+) l$- — 4-). 


(C.13) again shows that our propagating functional is properly normalized. 


Similarly, we expect the reduced density matrix to reduce to (C.1) when t > 0. Using 
(C.12) and a;;(t + 0) — 0, the exponent in (C.2) is 


Leg i l aa 1 Lig 5.4 119 call 
z A i i (C.14) 
OS 2 2 
ETEY 
where we used C + at. And the normalization factor in (C.2) takes the form 
aei i 
o E = —-. C.15 
ya a 4 &@ Qne OyT ( ) 


Combining (C.14) and (C.15), we reproduce (C.1). 
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